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\ Abstract 

Late-time entropy production from scalar decay arises in scenarios like thermal inflation, proposed 
to dilute long-lived, massive fields like the gravitino and the moduli. The scalar decay may continue 
into Mev-scale temperatures and affect BB nucleosynthesis. The effect of such entropy production 
on electron neutrino decoupling is studied. A lower bound of about 10~ 22 Gev is estimated for the 
scalar decay constant, such that, for higher values of the decay constant, standard electron neutrino 
^ ■ decoupling is unaffected. 
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1 Introduction 



Supersymmetric and string theory models throw up long-lived, massive fields like the grav- 
itino, the Polonyi, the moduli and the dilaton [1]. Stability of the corresponding particle 
may lead to over-abundance, while decay during baryogenesis or BBN nucleosynthesis may 
disturb 77 or nuclear abundances too much [2]. Thermal inflation [3] has been proposed as 
a method of diluting away troublesome fields. A scalar field, the flaton, is used to start 
inflation at a temperature of about 10 7 Gev. The inflation stops at a temperature of the 
order of the flaton mass, typically < 10 3 Gev. Such a scalar field goes on decaying into 
Mev-scale temperatures with potential trouble for BBN nucleosynthesis. 

The effect of flaton decay on the effective number of neutrino species, N e ff, and the 
neutrino distribution function during nucleosynthesis has been recently studied in detail 
[4, 5]. Here, we study the effect of scalar decay on decoupling of the usual left-handed, 
massless neutrinos. As the phenomenological details of such a decaying particle are still quite 
open, we consider the general effect of late time entropy production on electron neutrino 
decoupling, assuming that the decaying scalar once dominated the energy density of the 
universe. The electron neutrino decoupling temperature is a sensitive input of BBN, and 
we estimate lower bounds on the scalar decay constant, demanding that this decoupling 
temperature should not be disturbed. 

2 Entropy Generation and Decoupling 

The entropy generation rate is calculated from 

dS* = 

where is the scalar energy density at radiation temperature T, all massless particles 
being collectively called radiation. Let V be the scalar decay constant. Writing <E> = a 3 /9^, 
the equation for the evolution of the scalar energy density is = — r<3?, with the solution 
= $>Ee~ T ( t ~ tE ^[6], where <&e is the value of $ at a fiducial time t = tE- This gives 

5$ = {T/T)^ E e~ r(t - tE \ (1) 

What is the standard situation in the absence of scalar decay? In the era of neutrino 
decoupling, the electron neutrinos interact via charged and neutral currents, while muon 
and taon neutrinos interact via neutral currents alone. So, the muon and taon neutrinos 
decouple at a slightly higher temperature than the electron neutrinos. Then, the change in 
electron neutrino {y ) number occurs due to the process v + v — > e~ + e + . 

The decoupling of electron neutrinos is governed by the integrated Boltzmann equation 
[7, 8] 

h + 3Hn = - < a\v\ > (n 2 - n EQ ), (2) 

where n is the number density of the electron neutrinos, ueq their equilibrium number 
density, and < a\v\ > the thermally averaged cross-section times relative velocity. Assuming 
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an absence of neutrino degeneracy, we use the Boltzmann distribution function to give 

n EQ = T 3 /tt 2 . (3) 

Following the usual method of calculation [8, 9, 10], neglecting the masses of the neutrino 
and the electron, 

< a \v\ >= -G 2 F [(C Ve + I) 2 + (C Ae + lf]T 2 

4.112 x 10^° 2 

The decoupling may be taken to start when the following relation just holds. 

- < a\v\ > (n 2 - n 2 EQ ) < 3Hn. (5) 

Denoting (a 3 ) times entropy density as S, the net contribution, into the v, v sector, of 
the covariant divergence (l/a 3 )£ of the usual entropy density current, due to the process 
e~ + e + — ► v + u, is —2a < cr|?;| > (n 2 — n^g^ll], where 3 

n = n EQ e~ a . (6) 

If there is additional entropy generation due to scalar decay, there will be a corresponding 
extra term on the RHS of (2). To find a lower bound to T, such that scalar decay does 
not affect neutrino decoupling, we require this additional term to be small compared to the 
term — < a\v\ > (n 2 — n 2 E Q) on the RHS, at decoupling. To find the lower bound on T, we 
consider the worst case scenario when the entire additional entropy due to scalar decay is 
transferred 4 to the neutrinos. Then, if decoupling is to be kept undisturbed, we can require 
that the contribution, to the v, v sector, of the covariant divergence of the entropy current 
due to scalar decay should be less than that due to the usual annihilation process, or 

-^S$ < —2a < a\v\ > (n 2 — n 2 FO ), 
or v 

when (5) holds. This means that the required lower bound to T is to be obtained from the 
criterion 

-^5$ < 2a(3Hn), (7) 
a d 

3 In [11], the distribution function of the decoupling particle is taken as /(p) = e ~ a{t) ~ f3{t)E{v) (1 + 
For sufficiently strong elastic collisions, it is shown that and the process is well represented by the 

time-dependent parameter a(t). Despite the resemblance, a is not a chemical potential. It parametrises the 
distribution and is the same for particle and anti-particle, instead of being different in sign. 

The 2 factor in 2a arises because our equation (2) considers the change in particle density alone, while S 
considers net entropy transferred to both particle and antiparticle. [11] does not have the 2 factor as particle 
and anti-particle are there considered together as a single species. 

4 In current models, it is usually assumed that the scalar cannot decay directly to neutrinos [4]. Any 
massive decay products of the scalar quickly decouple and the entropy due to scalar decay passes to the 
e~, e + , 7 sector, and thence to the neutrinos via the annihilation process e" +e + — » v + v. As electromagnetic 
interactions keep the e~ ,e + ,y sector very near equilibrium, almost all the entropy may be supposed to pass 
to the neutrinos. 
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where a corresponds 5 to its decoupling value, in the absence of scalar decay, worked out 
from the equality corresponding to (5). 



3 Estimating the Entropy Generation 

(1) indicates that S§ depends on <&e- In the absence of definite phenomenological values, 
§E is to be estimated indirectly. We are interested in the era when radiation domination 
has begun, and (1) is dominated by the exponential. Then, the pre-exponential may be 
estimated approximately, to an order of magnitude. So, in the pre-exponential, we define 
the initial epoch £g by putting p^E = Pre-, Pr being the radiation energy density. Radiation 
domination is taken to start after this epoch. Further, we neglect entropy generation while 
estimating the pre-exponential in (1), and put 

§E = a\ P4>E = a% PRE = (vr 2 /30)^a|T| 

= (tt 2 / 30) g*a 3 T 3 T E , (8) 

where g*,a,T refer to the epoch of decoupling. 



3.1 Estimate of T E 

We assume that when the neutrinos decouple, p^ <C pr, such that p^/pr cannot be ne- 
glected, but its higher powers can. It is possible to show [12], in fact, that decoupling is 
not possible when p^ Pr- In the era p^ <C Pr of incomplete radiation domination, the 
Friedmann equation is put in the form 

3M 2 pi a 4y R h V ' 

where R = a A pR. If, well into this epoch, the correction term &a/R on the RHS of (9) is 
neglected, the full radiation domination relations are found: 

H = — , and 

2t' 

a = At^, (10) 

A being a constant. 

The <3? evolution equation has, as solution, a falling exponential in t, viz. $ ~ e~ rt . 
Instead of taking the falling exponential in t directly, a suitable approximation to the 



5 We assume that the electron neutrino distribution function continues to be parametrised as f(p) = 
e -a(t)-/3(t)E(p) even w ] leQ there is scalar decay. It is shown in meticulous detail in ref.[4] that, for low 
values of the scalar decay constant, the distribution is not thermalised. There is a deficit from the thermal 
distribution. So, our assumption is tantamount to an approximation of the deficit by a uniform factor e~ a , 
with a > 0. 
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correction term Qa/R on the RHS of (9) is first worked out. Let to be a sufficiently late 
epoch, when $ = <l?o ~ 0. Then, for use only in the correction term &a/R, one takes 

*-*° = *<7>-*^> = fl'°<K>- 

Neglecting <I>o,l/io compared to 3>, respectively, an approximation 

$«p (11) 

will be used only in the correction term <&a/R, i.e., in the correction term, the falling 
exponential will be approximated by a rectangular hyperbola. B is a constant. 
A similar approximation is considered for R. From 

d d 
-^[a 3 ( P(l) + p R )}+ PR —a 3 = 0, 

one obtains 

R = aT<S>. (12) 

It ought to be mentioned that R refers to the total radiation present, and not only to 
that produced by decay. However, the change in R is due to <fi decay and consequent 
entropy production. In the absence of this decay, R = 0. Using (10) and (11) in (12), and, 
integrating, one obtains approximately, for use only in the correction term <3?a/i?, 



R-R E ^ 2ABT(ti -t 2 



If ts is sufficiently early compared to t, and there is sufficiently copious radiation production 
since £e, it is sufficient to take 

R « 2ABTV2 

in the correction term <&a/R. (10) and (11) are now used to give, in the correction term, 

r 



X = H 



- la ^ 
Introducing the variable x in (9), we get (9) in the form 

H = T - 

X 

4.461 x 10~ 19 m2/ l,i 
= Gev ~x (14) 

putting g* = 10.75. (14) is a good equation for x S> 1. But, we approximate Tg in the 
pre-exponential term <E>£ of (1), by putting x = xe = 1 in (14) (corresponding to = pR, 
or $a = R, at t = *b) to get 

T E = i = (15) 

<JV2 x 4.461 x 10~ 19 /Gev 
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T D /Mev 


a 


x D 


T/(W-' M Gev) 


1 


3.4705 


17.19 


7.9 


2 


1.447 


14.36 


26.5 


3 


0.5644 


14.505 


60.2 



Table I: Neutrino a values at decoupling for different decoupling temperatures 
and the lower bound on the decay constant T 

4 Numerical Results and Conclusions 

Now, BB nucleosynthesis is consistent with a neutrino decoupling temperature Tp ~ a few 
Mev. Here, we take Tjj = 1, 2, 3 Mev, and find corresponding values of a, in the absence of 
scalar decay, from the equality corresponding to (5), using (3), (6) and H = 4.461 x 10 _19 T 2 , 
the last being the usual radiation domination value of H. The results are shown in the first 
two columns of Table I. 

Next, we consider scalar decay. The value of Tg from (15) is used in (8) and the latter is 
put into (1). Utilising (14) to relate the temperature and x where necessary, and estimating 
tE from H = l/(2i), (13), andx^ = 1, the equality corresponding to the relation (7) is now 
solved for x = xo, keeping the same values of a already found to correspond to decoupling 
temperatures Tjj = 1,2,3 Mev in the absence of scalar decay. The results are in the third 
column of Table I. Finally, V is found corresponding to these values of xd from (14) for 
Tjj = 1, 2, 3 Mev. The form of S$ and the inequality (7) show that these values correspond 
to a lower bound on T. 

We conclude that if the entropy produced by scalar decay is not to disturb the neutrino 
decoupling temperature, the lower bound on the scalar decay constant must be of the order 
of 8 x W~ 24 Gev , 26.5 x 10~ 24 Gev , 60 x W~ 24 Gev, corresponding to neutrino decoupling 
temperatures of 1,2,3 Mev. So, neutrino decoupling is unaffected, and BB nucleosynthesis 
safe on this count, if the scalar decay constant is greater than about 10~ 22 Gev, correspond- 
ing to a reheating temperature of 8.6 Mev. 
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